AO-A087  398 

UNCLASSIFIED 

!■<  I 


MASSACHUSETTS  INST  Of  TECH  CAMBRIDGE  LAB  FOR  INFORMA— CTC  F/6  10/2 
ROBUST  CONTROL  OF  MULTITERMINAL  DC/AC  SYSTEMS* (U) 

AU«  79  N  A  LEHTOMAKI*  S  M  CHAN*  N  R  SANOELL  DOE-ET-70-C-O 1-3395 


8  6  €  1 8  o  v  av  j^joo  m  3QQ 


rr. 


3j3  \  UM 

r 


7  L/ 


ROBUST  CONTROL  OF  MULTITERMINAL  DC/AC  SYSTEMS# 
S-  —  —  ' 

by 


N.A.  Lehtooaki ,  S.M.  Chan,  N.R^^andell.Jr. /H.^Athans,  J.  fcarrig/ 
Laboratory  for  Information  and  Decision  Systems 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

C.E.  Grund,  F.  Nozari,  R.  Hauth 
— s —  .  Rlactrlc  Utility  Systems  Engineering  Department 

t  7 /  — L  9^7  General  Electric  Company 


—  / 


Schenectady,  New  York  12345 


INTRODUCTION 


In  the  past  two  decades,  dc  power  transmission  has 
provided  a  viable  alternative  to  ac  in  special  applica¬ 
tions  such  as  long-distance,  point-to-point  transmis¬ 
sion,  asynchronous  connection  between  power  systems  of 
different  frequencies,  and  interconnection  of  systems 
via  underground  and  undersea  cables.  Common  to  all 
these  applications  is  the  two-terminal  dc  link  which 
consists  of  a  rectifier  and  an  inverter  connected  by  a 
dc  transmission  line.  The  main  advantage  of  dc  trans¬ 
mission  lies  in  its  ability  to  transmit  a  large  block 
of  power,  with  very  little  los3,  from  one  part  of  the 
network  to  another  without  the  inherent  synchronous 
and  stability  limitations  of  ac  systems .  Although 
most  dc  Units  to  date  are  built  for  the  purpose  of 
bulk-power  transmission,  it  was  discovered  recently 
that  by  modulating  the  dc  current,  a  dc  link  can  be 
used  as  an  extremely  effective  control  element  for  damp¬ 
ing  some  of  the  inter- area,  electromechanical  oscil¬ 
lations  in  an  integrated  ac/dc  system. 

— ==r-~i.The  successful  application  of  two  terminal  dc 
links  suggests  that  even  greater  flexibility  in  power 
dispatch  and  enhancement  in  stability  can  be  realized 
by  a  multitenninal  dc  (MTDC)  system  embedded  in  an 
Interconnected  ac  grid.  Several  utilities  are  looking 
into  the  possibility  of  building  MTDC  systems  by  adding 
new  terminals  to  the  existing  dc  links  or  by  connecting 
dc  links  that  are  located  in  close  proximity.  It  is 
quite  conceivable  that  MTDC  will  become  a  major  part 
of  our  energy  delivery  system  by  the  end  of  this 
century .  - 

A  parallel  connected  MTDC  system  is  envisaged  as 
a  network  of  n  (n>2)  inverters  and  rectifiers  intercon¬ 
nected  by  dc  transmission  lines.  In  the  steady  state, 
n-1  converters  are  operated  under  a  constant-current 
mode,  and  the  nth  converter  with  the  lowest  ac  com¬ 
mutating  voltage  (voltage  that  is  applied  to  the  ac 
side  of  a  converter)  is  operated  under  constant  delay- 
angle  or  constant  extinction-angle  mode  and  determines 
the  voltage  of  the  entire  MTDC  system.  In  a  transient 
condition  the  terminals  of  an  MTDC  system  can  poten¬ 
tially  be  modulated  in  a  coordinated  manner  to  damp,  out 
the  inter-area  electromechanical  oscillations. 

The  design  of  the  MTDC  control  system  offers  a 
unique  challenge;  since  it  involves  a  large-scale, 
multi variable  system  where  the  sensors  and  actuators  are 
situated  far  apart,  and  the  parameters  of  the  system 
are  difficult  to  estimate  and  may  change  abruptly  due 
to  failures  such  as  short  circuits  or  generators  out¬ 
ages.  A  workable  design  for  such  a  system  must  be 
sufficiently  robust  to  tolerate  actuator  and  sensor 
failures,  urunodelled  dynamics,  and  changes  in  system 
parameters.  Moreover,  economic  and  reliability  consid¬ 
erations  dictate  that  some  form  of  decentralized 
feedback  scheme  be  used,  if  possible. 
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'  ^  The  intent  of  this  paper  is  to  illustrate  how  the 

MTDC  control  problem  can  be  analyzed  from  a  physical 
/  viewpoint,  and  how  the  physical  insights  can  be  uti¬ 
lized  to  arrive  at  a  reasonable,  robust  control 
design.  The  main  emphasis  is  on  robustness.  Some  of 
the  latest  results  in  robustness  theory  are  used  in 
the  design  stage  as  well  as  in  analyzing  the  final 
design.  £. _ 

Although  this  paper  presents  a  general  design 
methodology  for  an  MTDC  system,  a  hypothetical,  five- 
terminal  dc/ac  system  is  used  as  an  example  to  illus¬ 
trate  the  steps  of  the  design  process.  The  model  of 
this  system  is  described  in  Section  I, -and  a  physical 
interpretation  of  the  design  model  from  a  modal  pers¬ 
pective  is  given  in  Section  II.  Section  III  shows 
hew  the  physical  insights  can  be  used  to  construct  a 
state  penalty  matrix  for  an  LQ  design,  and  the  concept 
of  gain  observability  is  utilized  to  study  the  trade¬ 
off  between  the  LQ  centralized  design  and  other  decen¬ 
tralized  schemes.  The  robustness  of  the  control  design 
is  studied  in  Section  TV,  after  a  brief  review  of  a 
relevant  theorem  in  robustness  theory.  Finally,  some 
results  of  the  open-  and  closed-loop  systems  are 
compared  in  Section  V. 

I.  MODELLING 

The  power  system  being  studied  consists  of  an  n- 
terminal  MTDC  system  embedded  in  an  ac  system  with  m 
groups  of  coherent  machines,  each  of  which  is  repre¬ 
sented  by  the  classical  constant-voltage-behind- 
transient-reactance  generator  model.  For  each  gener¬ 
ator  the  electrical  angle  6  is  governed  by  Newton's 
Second  Law  for  a  rotating  mass : 

jS  +  Dfi  -  P  -  P„  (1) 

nx  q 


where 

J  is  a  constant  related  to  the  inertia  of 
the  generator, 

p  is  the  mechanical  power  supplied  by  the  prime 
mover, 

P  is  the  electrical  power  injected  into  the 
power  network,  and 

D  is  a  constant  related  to  the  damping  of  the 
machine. 

It  is  further,  assumed  that  the  mechanical  power  re-  ' 
mains  constant  over  the  time  period  of  interest. 

The  dynamics  of  the  dc  converter  and  the  dc  net¬ 
work  arc  ignored,  because  they  are  much  faster  than 
the  electromechanical  dynamics  being  studied.  If  the 
dc  terminals  are  viewed 
that  the  control  action 
all  time. 
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as  actuators,  this  implies  /  j 
is  equal  to  the  command  at 


i  ijnoy' 
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Under  these  assumptions,  a  linearized  model  is 
obtained  by  first  solving  for  the  steady-state  solu¬ 
tion  of  the  ac/dc  system,  and  then  linearizing  the 
equations  about  the  steady-state  point.  The  result  is 
a  linear  time-invariant  model  of  the  form 

x(t)  -  hjc (t)  +  Bu(t)  (2) 

where 

A  T 

x  -  (A6, ,  Aw, , . . .  ,A5  ,Aw  )  are  the  changes  in 

—  x  l  mm 

electrical  angle  and  frequency  of  the  m 
generators , 

A  T 

u=  (Alj,...,  Al^)  are  the  changes  in  current 
injected  by  the  dc  terminals . 

(The  current  injection  of  the  first  dc  terminal  is 
not  a  control  variable) .  A  and  B  are  matrices  of 
appropriate  dimension. 

In  the  subsequent  discussion,  the  "A”  symbol  will 
be  dropped,  and  all  quantities  are  deviations  from  the 
steady-state  solution  unless  otherwise  specified. 

It  is  important  to  emphasize  that  the  design  model  (2) 
does  not  take  into  account  the  dynamics  associated 
with  the  excitation  system,  power-system  stabilizer, 
shaft  torsional  vibrations,  and  a  host  of  other  phjs- 
ical  phenomena  which  are  significant  at  frequencies 
greater  than  30  rad/sec;  thus  the  design  model  (2)  is 
valid  only  for  low  frequencies.  This  fact  will  be 
Important  when  designing  the  feedback  system. 

For  the  purpose  of  illustration,  a  hypothetical 
ac/dc  system  shown  in  Figure  1  is  used  throughout. 

A  physical  interpretation  of  the  dynamics  associated 
with  this  system  will  be  given  after  a  brief  discus¬ 
sion  on  the  methods  of  modal  analysis. 

II.  CONTROLLABILITY,  OBSERVABILITY  AND  PHYSICAL 

INTERPRETATION  OF  MOOES 

Since  the  linearized  model  on  which  the  control 
design  is  based  has  limitations,  as  all  models  do,  it 
is  important  to  be  able  to  physically  interpret  the 
system  modes  in  order  that  the  control  system  design 
may  be  evaluated  in  part  by  how  it  affects  them.  It 
fs  also  Important  to  assess  beforehand  the  inherent 
ability  of  the  available  control  inputs  to  accomplish 
the  performance  objective  using  only  information  fur¬ 
nished  by  the  available  output  measurements.  This  is 
done  by  computing  the  open-loop  system  eigenvalues 
the  corresponding  left  and  right  eigenvectors,  deter¬ 
mining  what  each  mode  physically  represents  and  then 
determining  the  controllability  and  observability  of 
each  mode. 


The  reason  for  the  adjective  "right"  and  "left"  ia 
obvious  by  noting  that 

AV  -  VA  (6 

WA  -  AW  O 


A  =  diag(X, ,X  , ...,A  1  . 

l  *  n 


„  (8) 


Using  the  right  and  left  eigenvectors,  the  matrices 
At 

A  and  e  may  be  written  as  dyadic  expansions: 


A  “  l  <9> 

i=l  1 

At  ?  V  T 

e  -  l  e  1  (10) 

i*>l 

The  unforced  response  of  the  system  is  new  given  by 
n  At 

x(t)  -  l  (e  1  v, )  (wT  x)  (11) 

i=l 

which  shows  that  x(t)  is  simply  a  linear  combination 

of  n  tine-varying  vectors,  e  v^,  or  modes  each 
weighted  by  the  scalar  w^  Xg.  A  physical  interpre¬ 
tation  of  a  right  eigenvector,  v^,  can  best  be  illus¬ 
trated  graphically  on  an  Argand  diagram  where  the 
components  of  v.  are  plotted  as  vectors  in  the  complex 

V 

plane.  The  quantity  e  v^  can  be  visualized  as 
rotating  all  components  of  at  the  same  rate 
ImlX^  and  simultaneously  shrinking  or  expanding  them  at 
the  rate  exp {Re [X^t] }  depending  on  whether  X^  is  sta¬ 
ble  or  unstable.  This  is  similar  to  th#  phasor  re¬ 
presentation  of  sinusoidal  signals  and  is  depicted  in 
Figure  2 . 

The  left  eigenvectors,  on  the  other  hand,  deter¬ 
mine  a  linear  combination  of  states  that  oscillate  at 
the  complex  frequency  equal  to  the  corresponding 
eigenvalue  since 


§£■  {  w^x(t)}  -  wj  Ax(t)  »  \l  {wj  x(t)}  . 


Modal  Analysis 


Physical  Interpretation  of  System  Modes 


To  facilitate  the  subsequent  discussion  on  the 
physical  interpretation  of  the  system  modes  the  role 
of  the  left  and  right  eigenvectors  of  the  system 
matrix  A  is  reviewed.  Consider  the  open-loop  system 

x(t)  =  Ax(t)  i  x(0)  “  x^  (3) 

where  the  eigenvalues  of  A,  X^ .X^ , . . . .X^ ,  are  assumed 

to  be  distinct.  Let  v,  and  wT  be  the  riqht  and 
— 1  —1 

left  eigenvectors  respectively  of  A  associated  with 
the  eigenvalue  X^^  and  define  the  matrices 


The  open-loop  eigenvalues  of  the  system  are 
plotted  in  Figure  3.  The  Argand  diagram  associated 
with  each  of  the  four  pairs  of  complex  eigenvalues 
reveals  that  the  frequencies  of  certain  areas  are  ap¬ 
proximately  180*  out  of  phase  with  that  of  the  other 
areas.  This  is  interpreted  as  an  electromechanical 
oscillation  in  which  two  coherent  groups  of  areas  for 
a  given  mode  swing  against  each  other.  For  example, 
in  Figure  4  for  the  -,094+j3.80  modes,  the  areas  i  and  2 
swing  against  areas  3,4  and  5  with  areas  2  and  S 
being  dominant  in  the  oscillation.  The  two  remain¬ 
ing  eigenvalues  are  real.  For  the  eigenvalue  at  zero 
the  right  eigenvector  is 


Ax  -  A(x  +  6^)  .  (14 

In  other  words,  only  the  difference  in  angle  states 
is  important  since  any  constant  may  be  added  to  all 
the  angle  states  without  affecting  the  dynamic  res¬ 
ponse  of  the  system.  The  left  eigenvector  for  the 
remaining  mode  at  -.18  has  the  form 

v£  -  (0,H1(0,H2,0,H3,0,H4,0,H51  (15 

where  is  the  intertia  of  area  1  divided  by  the 
total  system  intertia.  Equation  (12)  implies  that 

u(t)  -  -.18w(t)  (16 


u(t)  =  wT  x(t) 


Z  V 

i»l 


is  interpreted  as  the  average  frequency  of  the  power 
system. 

The  foregoing  analysis,  clearly  demonstrates  the 
role  of  the  left  and  right  eigenvectors  in  determining 
the  physical  interpretation  which  may  not  be  immediate¬ 
ly  apparent  for  some  of  the  system  eigenvalues . 

The  left  and  right  eigenvectors  are  also  very  useful 
in  determining  the  controllability  and  observability 
of  the  modes. 

Modal  Controllability  and  Observability 


puts  and  n  states 

x(t)  v  ax (t)  +  Butt),  x(0)  *=  x 

—  ■  ""V 

£(t)  -  Cx(t) 

Using  the  variation  of  constants  formula  gives 


...  At 
x(t)  »  e  Xjj 


_y  (t)  -  CeAtxQ 


>  ye 

0 

t 

'*/ 


Bu(T)dT 


CeA<t'T)Bu(T)dT 


For  linear  operators  L(-)  and  H(')  let  R  (L)  denote 
the  range  of  L(-)  and  N(H)  denote  the  nullspace  of 
H(«).  The  system  defined  in  (18)  and  (19)  is  com¬ 
pletely  controllable  if  and  only  if 


/ 


Bu(l)dT 


and  is  completely  observable  if  and  only  if 
W(H)  -  {0} 


if  the  state  of  the  system  can  be  driven  anywhere  in 

Rn  by  some  control  u(T)  applied  from  T-0  to  t. 
Analogously,  (24)  implies  that  given  the  output  ^(t) 
and  the  input  u(T)  for  0<r£t,  the  initial  state,  and 
thus  the  state  trajectory,  can  be  determined  or 
observed. 

There  are  different  ways  of  determining  whether 
the  conditions  (22)  and  (24)  are  satisfied.  However, 
the  definitions  of  controllability  and  observability 
specified  in  these  conditions  are  qualitative  con¬ 
cepts  and  do  not  provide  quantitative  measures  of  how 
controllable  or  how  observable  a  system  may  be;  they 
can  only  say  what  part  of  the  system  is  controllable 
or  observable  and  what  part  is  not.  Part  of  the  prob¬ 
lem  is  that  it  is  not  immediately  apparent  what  is 
meant  by  terms  like  "strongly  controllable"  or 
"almost  observable",  since  it  is  possible  to  define 
them  differently  depending  on  their  intended  use. 

One  possible  way  to  quantify  controllability  and 
observability  is  to  measure  the  change  with  respect 
to  perturbations  in  the  nominal  system  model.  In 
other  words,  if  "small"  perturbations  in  the  nominal 
system  change,  the  controllability  of  certain  states, 
those  states  would  be  considered  "weakly  controllable* 
or  "almost  uncontrollable." 

A  different  but  not  unrelated  way  of  measuring 
the  controllability  and  observability  is  based  on  the 
intuitive  idea  of  how  strongly  the  modes  can  be 
excited  by  the  inputs  and  how  strongly  the  modes  ap¬ 
pear  at  the  outputs,  respectively.  This  is  made 
more  precise  by  considering  the  original  qualitative 
definitions  of  controllability  and  observability  in 
terms  of  the  modes  of  the  system.  Consider  L(u)  and 


p  out- 

H(Xg)  in  terms  < 

(18) 

n  1 

(19) 

L(u)  *=  l 
i-1 

n 

HIXq)  "  l 

1M  • 


From  (26)  it  is  clear  that  if  W ,  B-0  for  scow  i, 

x 

then  L(u)  is  a  linear  combination  of  at  most  n-1 
eigenvectors  for  any  u>  hence  dim  R (L) <  n  and  the  sys¬ 
tem  is  not  completely  controllable.  Similarly  from 
(27),  if  Cv^-0  for  some  i,  then  HCv^J-O  but  v^O; 

hence  dim  N (H) >  0  and  the  system  is  not  cotapletely 
observable.  These  facts  are  immediately  obvious 
once  ^(t)  is  written  as 

r  f  T  T  f  "V  1  M 

£(t)  «  I  Cv^  w^b  J  e  uCtMtJ®  1  .  (28) 

T  th 

If  w.B*0  then  the  input  u  cannot  excite  the  i- — mode 

and  if  Cv.=0  the  i — mode  does  not  appear  in  the 

*  th 

output  y_(t) .  In  these  cases  th#  i —  mode  is  said  to 

be  uncontrollable  or  unobservable  respectively.  If 
the  matrices  B  and  C  are  given  by 

B  -  lb  ,b  . b  1  (29) 

1  t  Ul 


.  ,  ,  6  _  At 

HUq)  -  Ce  ^ 


c  “ 


th  T 

where  b.  is  the  i —  column  of  B  and  c  is  the  i- 
— i  — i 

row  of  C  then  the  matrices  WB  and  CV  given  by 


Tv 

T 

T, 

-1-1 

—1—2  ‘ - • 

.  -l-m 

..Tv 

wTb  . . . . 

T, 

-2-1 

-2-2 

.  — 2— m 

bars  of  unit  height.  In  Figure  6  the  bar  heights  have 
been  scaled  coluran-wise  such  that  for  a  given  mode  the 
bars  of  the  input  and  output  from  which  the  mode  is 
most  controllable  and  observable  have  unit  height. 
Notice  that  the  arbitrary  phase  reference  mode  is  not 
observable  from  any  frequency  output  but  is  control¬ 
lable  from  every  current  input.  From  both  Figures  5 
and  6  it  is  evident  that  the  most  observable  oscil¬ 
latory  modes  in  each  frequency  output  are  also  the 
most  controllable  from  the  input  of  the  same  area. 

This  has  important  ramifications  for  the  control  sys¬ 
tem  design  in  that  the  sensors  and  actuators  may  be 
colocated,  obviating  the  need  for  expensive  ~onmunica- 
tion  links. 


have  a  nice  interpretation:  the  magnitude  of  the 

T  th 

entry  w.b  of  WB  measures  how  much  the  j —  input 

th 

of  u  affects  the  i —  mode  and  the  magnitude  of  the 

T  th 

entry  c, v,  of  CV  measures  how  much  the  j —  mode  ap- 

~i~i  th 

pears  in  the  i —  output  of  (t) . 

If  the  rows  of  C  and  the  columns  of  V  are  normal¬ 
ised  to  have  unit  euclidean  length ,  the  magnitude  of 
T 

£.v^  represents  the  cosine  of  the  angle  between  the  vec 

tors  c  and  v.  .  Thus  when  c,  and  v  are  nearly 

‘  i  K  ^  ^  ^ 

orthogonal  the  k —  mode  Is  unobservable  from  the  i — 
output  for  some  small  perturbations  of  the  nominal 
system. 

In  interpreting  (31)  and  (32),  it  is  important 

to  keep  in  mind  that  the  magnitude  of  |  |wb  |  |  and 

•  X  *  ^ 

|  |c^v_j  |  1  depends  on  the  choice  of  units  of  the  input 

and  output  as  well  as  the  scaling  of  the  eigenvectors . 
If  the  units  of  the  input  and  output  are  considered 


fixed,  the  quantity 
making  ||wTb||  sma 


can  only  be  made  large  by 
“1  T 

t  V=W  and  w ,  v ,  *»  6 .  ,  . 


-IT 

small,  since  V=W  and  w.v, 
-1-j 


ftius  it  seems  reasonable  to  balance  the  left  and  right 
eigenvectors  corresponding  to  the  same  eigenvalue  by 
making  their  norms  equal.  This  ensures  that  a  mode 
does  not  appear  very  observable  but  not  very  control¬ 
lable  or  vice  versa  only  because  of  the  scaling 
chosen.  It  is  also  worth  noting  that  the  matrix  prod¬ 
ucts  CV  and  WB  are  invariant  under  similarity  trans¬ 
formations  of  the  original  system,  and  the  terms 

Cv  vTx  and  Cv.v?B  in  (28)  are  independent  of  the 
— i~i~0  — i-r 

scaling  of  the  right  and  left  eigenvectors  as  long  as 

>j*:V  For  the  power  system  model,  the  matrices  CV 
and  WB  (where  the  left  and  right  eigenvectors  are  bal¬ 
anced)  ,  have  the  magnitudes  of  their  elements  plot¬ 
ted  on  bar  graphs  in  Figures  5  and  6.  The  bars  are 
plotted  in  pairs,  the  one  on  the  left  for  the  obser¬ 
vability  at  the  indicated  output  and  the  one  on  the 
right  for  the  controllability  at  the  indicated  input. 
This  is  done  for  each  mode  of  the  system.  In  Figure 
5  the  heights  of  the  bars  in  a  given  row  have  been 
scaled  such  that  the  most  controllable  and  most 
.  V  .  . .vi  ~  .Vo  o ;  irm  i  nnnt  and  output  have 


3  AD  RATIO  REGULATOR  DESIGN 


The  basic  design  methodology  adopted  here  la  that 
of  using  a  linear  quadratic  regulator  to  obtain  a  first 
iteration  centralized  design  that  moves  the  open  loop 
poles  of  the  system  to  an  appropriate  closed  loop  pole 
region.  This  region  is  determined  approximately  on 
engineering  judgement  by  how  large  the  closed  loop 
bandwidth  may  be  made  without  allowing  unmodeled  high 
frequency  disturbances  or  dynamics  to  destabilize  the 
system.  The  centralized  design  is  done  with  an  eye 
towards  later  decentralization  which  uses  only  the 
measurable  frequency  states.  The  decentralized  control 
law  then  attempts  to  mimic  the  centralized  control  law 
by  including  only  the  important  feedback  gains  from 
the  centralized  feedback  gain  matrix. 

Centralized  Design 

The  optimal  feedback  control  law  that  minimizes 
the  cost  functional 


J(u)  «  /xr(t,Cx(t,  +  uT(t)Ru(t)i  Q>0)  R>0  (33) 


is  given  by  the  familiar 
u»(t)  ■=  -R_1BTKx(t) 


ATK+KA-VQ  -  KBR_1£k  -  0  (35) 

1/2 

under  the  usual  assumptions  that  (A;Q  )  is  detect¬ 
able  and  (A,B)  is  stabilizable.  The  state  weighting 
matrix  Q  is  selected  so  that  neither  the  real  average 
frequency  pole  nor  the  arbitrary  phase  reference  pole 
at  zero  are  moved.  It  is  not  desirable  to  move  the 
pole  at  zero  since  it  is  only  observable  from  the 
angle  states  which  are  not  measurable.  The  average 
frequency  pole  is  not  to  be  moved  since  the  load 
frequency  control  loop  is  to  control  that  mode.  The 
only  modes  that  are  to  be  changed  are  the  inter-area 
oscillatory  modes  since  they  require  significant  damp- 
*n9-  There  are  at  least  two  ways  that  the  state 
weighting  matrix  Q  may  be  selected  to  do  this.  The 
first  method,  a  modal  one,  is  more  general  than  the 
second  method  which  depends  on  the  physical  inter¬ 
pretation  of  the  system. 

Modal  Selection  of  0  Matrix 


Consider  the  diagonalization  of  the  system  (18) 
by  defining  a  modal  vector  £(t)  given  by 

£(t)  =  Wx  (t)  (3« 

which  results  in 

z(t)  -  Az (t)  ♦  NBu(t)  O' 


Atits  Au«d  in  uie  seieCLAun  ul  y  la  lu  weiyug 

only  those  components  of  £(t)  in  the  cost  that  cor¬ 
respond  to  poles  of  the  open- loop  system  that  are  to 
be  moved.  The  rest  of  the  modes  of  the  system  are 
made  cost  unobservable.  Thus,  if  D  is  a  diagonal 
matrix  of  weightings  for  each  mode  then  (33)  becomes 


J(fi) 


Dz(t)||2  +  uT(t)Ru(t)dt 
2 


(38) 


Note  that  the  term  ||D£(t)j|2  can  be  converted 
T  2 

into  x  (t)Qx(t)  with  a  real  Q  matrix  in  the  following 

way.  Define  T  as 


r  -  dw  -  rR  t  ivx 


(39) 


then 

||Dz(t)!|2  »  xT(t)WBD2Wx(t)  =  JcT(t)Qx(t)  (40) 

where 


average  frequency  deviation  weighting  for  the  i — 
area.  This  Q  only  weights  frequency  variables  so 
that  the  zero  pole  will  not  be  moved  and  since  only 
deviations  from  the  average  frequency  are  weighted 
the  average  frequency  is  used  as  a  reference  and 
hence  the  pole  at  -.18  will  not  be  moved.  Using  this 
method  also  requires  the  stabilization  of  the  aero 
phase  reference  mode .  Therefore  in  the  centralized 
design,  ei  is  subtracted  from  the  system  X  matrix  and 
is  used  as  the  system  matrix  when  solving  the  Riccati 
equation  (35) a  the  resulting  feedback  control,  however, 
is  applied  to  the  original  unstabilized  system.  Both 
methods  of  selecting  Q  require  only  a  few  iterations 
of  selecting  modal  weightings  or  frequency  deviation 
weightings  to  obtain  the  approximate  damping  of  the 
oscillatory  modes  required.  It  is  interesting  that 
both  methods  give  approximately  the  same  feedback 
gains.  The  centralized  design  that  is  used  throughout 
the  remainder  of  the  paper  is  based  on  the  weighting 
of  deviations  from  average  frequency.  It  should  also 
be  mentioned  that  the  control  weighting  matrix  R  was 
selected  as  the  identity  matrix  multiplied  by  a  scalar 
to  tradeoff  state  and  control  energy. 


fR  rR  +  ^ 


(41) 


There  is  a  slight  difficulty  with  this  method, 
however,  in  that  (A,Q  '  ]  may  not  be  detectable. 

This  is  the  case  because  the  arbitrary  phase  reference 
mode  at  zero  is  made  cost  unobservable .  One  way 
around  this  problem  is  to  use  a  reduced  model  that 
does  not  have  this  mode.  An  alternative  way  is  to 
artificially  stabilize  the  cost  unobservable  unstable 
modes  leaving  the  other  system  modes  fixed.  This 
new  system  is  essentially  the  same  as  the  old  system 
for  the  purposes  of  design  of  a  feedback  controller 
since  the  controller  will  not  attempt  to  stabilize 
the  unstable  cost  unobservable  modes.  This  is  ac¬ 
complished  utilizing  the  dyadic  expansion  of  the  A 

matrix  in  (9) .  Let  A  be  the  artificially  stabilized 
system  matrix  which  is  defined  by 


-  I 

USt 


e  L 
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(42) 


where 


Feedback  Structures 


In  order  to  describe  several  controllers  that 
result  from  modifications  of  the  centralized  control¬ 
ler  it  is  necessary  to  classify  several  different 
feedback  structures.  They  are  the  follcwing: 

(1)  LQ  Centralized  -  full  state  feedback 

(2)  LQ  Full  State  without  Angles  -  frequency  feed¬ 
back  only 

(3)  LQ  Major  -  local  area  frequency  and  frequency 
of  area  1  available  for  feedback  (the  major 
centralized  gains) 

(4)  LQ  Decentralized  -  only  local  area  frequency 
available  for  feedback. 

In  each  case  the  feedback  gain  for  the  controller 
structure  is  obtained  by  simply  setting  the  gains  in 
the  centralized  gun  matrix,  that  correspond  to  states 
not  available  for  feedback,  to  zero.  The  closed-loop 
eigenvalues  for  cases  1, 3  and  4  are  given  in  Figure  7. 
To  assess  how  closely  each  of  the  feedback  controllers 
mimic  the  centralized  feedback  controller  the  concept 
of  gain  observability  is  introduced. 


0  “  {i:  mode  is  unstable  and  cost 

unobservable)  ,  (43) 

Now  the  matrix  A  has  the  same  eigenvectors  as  A 
and  the  same  eigenvalues  as  A  except  for  those  as¬ 
sociated  with  the  index  set  (5.  This  ensures  that  the 
old  and  new  system  will  behave  the  similarly  from  the 
standpoint  of  controller  design.  A  somewhat  simpler 
variation  of  this  approach,  which  works  if  €  is  suf¬ 
ficiently  small,  is  to  define  A  as  A=A-ei.  Then  A 
ha3  the  same  eigenvectors  as  A  and  eigenvalues  e  . 
away  from  those  of  A. 

Average  Frequency  Deviation  Selection  of  Q  Matrix 

The  second  method  of  selecting  the  state 
weighting  matrix  Q  is  to  weight  only  each  machine's 
or  area's  frequency  deviation  from  the  average  fre¬ 
quency.  Thus  the  state  weighting  is 

xT(t)Qx(t)  S  £  a  (iu,  — <o) 2  (44) 

i-1  1 


Gain  Observability 

The  ability  to  realize  a  centralized  feedback 
law  in  a  decentralized  manner  depends  critically  on 
the  nature  of  the  centralized  feedback  law.  It  may 
not  even  be  possible  if  the  centralized  feedback  law 
depends  substantially  on  states  not  available  for 
feedback.  Thus  only  modes  observable  from  the  fre¬ 
quency  outputs  are  moved  in  the  centralized  design. 
Eliminating  gains  in  the  centralized  gain  matrix  that 
do  not  significantly  affect  the  feedback  is  commonly 
done  by  setting  the  "small"  gains  to  zero.  What  is 
considered  "small"  must  be  considered  in  reference  to 
the  expected  size  of  the  variations  in  the  state. 
However,  it  is  not  always  clear  what  the  expected  vari¬ 
ations  are  since. in  general,  some  of  the  state 
variables  may  not  be  physical  quantities.  Also,  while 
a  particular  gain  may  not  by  itself  be  insignificant, 
a  group  of  certain  gains  may  collectively  have  very 
little  effect  on  the  feedback  control  signal  imple¬ 
mented  since  cancellation  of  the  contributions  of 
different  gain-state  products  may  occur  when  summed. 
This  later  idea  is  made  precise  by  the  notion  of  gain 
observability.  Gain  observability  merely  determines 
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the  modes  being  fedback  In  each  control  channel  by 
using  the  gain  matrix  as  an  output  measurement  matrix 
and  computing  the  matrix  GV  where  G  is  the  gain  matrix 
and  V,  given  by  (4) ,  is  the  matrix  of  right  eigenvec¬ 
tors.  By  computing  GV  for  different  gain  matrices  G, 
the  modes  being  fedback  in  each  channel  may  be  com¬ 
pared.  This  is  done  in  Figure  8  for  che  four  feed¬ 
back  structures  previously  described.  From  Figure  6 
it  is  apparent  that  the  LQ  Centralized  and  LQ  Full 
State  without  Angle  feedback  cases  are  practically 
identical  indicating  it  is  not  important  to  try  and 
measure  the  angle  states.  The  LQ  Major  case  is  also 
quite  similar  to  the  LQ  Centralized  case  while  the 
LQ  Decentralized  case  is  somewhat  less  so.  Thus  gain 
observability  indicates  the  dominant  modes  in  each 
feedback  channel  and  provides  insight  on  how  to  cons¬ 
truct  a  decentralized  or  limited  state  feedback  con¬ 
troller  that  approximates  a  centralized  controller. 

As  a  matter  of  practical  importance,  each  control 
signal  in  the  real  system  is  passed  through  a  washout 
which  washes  out  any  constant  frequency  errors  that 
may  be  generated.  The  washouts  have  a  transfer 

function  — ~  and  are  included  in  the  subsequent 
8+  .  2 

analysis  as  part  of  the  controllers. 

IV.  ROBUSTNESS 


In  any  control  system  design  there  is  a  degree  of 
uncertainty  abouv.  the  model  parameters  and  about  the 
nature  of  the  disturbances  that  act  as  exogenous  in¬ 
puts  to  the  model.  A  well  designed  control  system 
must  be  able  to  tolerate  such  uncertainty  without 
destabilizing  the  system.  Classically,  a  measure  of 
robustness  has  been  the  notions  of  gain  and  phase 
margin  for  SISO  feedback  systems.  In  the  multivari¬ 
able  version  of  the  Nyquist  stability  theorem,  the 
gain  and  phase  margins  do  not  provide  an  adequate 
measure  of  the  robustness  of  the  nominal  feedback  sys¬ 
tem  since  arbitrarily  small  simultaneous  perturbations 
in  the  nominal  system'roay  cause  instability  of  the 
closed- loop  system  even  though  there  are  good  gain  and 
phase  margins  in  each  feedback  channel. 

The  multivariable  Nyquist  criterion  counts  the 
encirclements  of  the  origin  by  the  det(I+C(j<i>) )  where 
G(s)  is  the  loop  transfer  functions.  In  order  to 
Change  the  nunber  of  encirclements  of  the  origin, 
det(I+G(ju))  must  be  zero  for  some  w.  Ulus  for  nom¬ 
inally  stable  closed- loop  systems  the  det(I+G)  trie3 
to  measure  the  "distance"  of  I+G  to  the  critical  point 
but  Is  not  a  good  measure  because  det(I+G)  may  be 
very  sensitive  to  small  changes  in  G.  In  the  scalar 
case  det(I+C)  becomes  1+g  and  small  changes  in  g 
yield  only  small  changes  in  the  Nyquist  diagram  of 
1+g  or  g  so  that  this  problem  does  not  arise.  There¬ 
fore,  in  the  multivariable  setting  what  is  required 
is  an  insensitive  measure  of  the  "distance"  of  I+G  to 
the  origin.  The  minimum  singular  value  of  I+G  is 
such  a  measure.  It  measures  the  minimum  size  of  the 
additive  perturbation,  AG,  required  to  make 
det (I+G+AG)=0  or  I+G+AG  exactly  singular.  The  min¬ 
imum  and  maximum  singular  values  are  denoted  as  0 
and  0  respectively  and  are  defined  by 


min 

x/0 

1  1*5112 

1 

,  if  A  ^  exists’ 
(45) 

l|x|l2 

1  |  1 1  1 2  ' 

max 

>1*5' >2  A 

IIa||2 

(46) 
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nominal  system  that  will  make  the  closed— loop  system 
unstable.  It  has  been  known  for  some  time  that  linear 
quadratic  regulators  for  single  input-single  output 
systems  have  inherently  good,  guaranteed  gain  and  phase 
margins,  of  namely  infinite  upward  and  50%  reduction 
gain  margins  and  +60*  phase  margin  tl] .  This  is 
because  for  LQ  regulators  1 1+g  ( jw)  |>1  Vw  so  that  the 
Nyquist  diagram  of  g  must  avoid  the  unit  disk  centered 
at  (-1,0) .  This  is  shown  in  Figure  9.  For  multi- 
Vuiable  systems,  Kalman’s  well-known  inequality  for 
LQ  regulators  [1,4]  generalizes  with  R*l  to 

(I+G(-s))T(I+G(s))>  I  Vs  (47) 

so  that 

0(I+G(jw))>  1  VU  (48) 

is  analogous  to  the  scalar  case  whefe  |l+g  (jw) )  |  >  1. 

It  has  more  recently  been  shown  [2]  that  for  diagonal 
R>0  and  Q>0  that  the  multivariable  LQ  regulators  have 
the  same  gain  and  phase  margins  in  each  feedback 
channel  simultaneously  as  in  tha  single  input-single 
output  case.  However,  as  is  well  known,  the  infinite 
gain  margin  is  not  a  reality  for  the  real  world  be¬ 
cause  the  nominal  model  is  only  valid  within  a  limited 
frequency  range  and  furthermore  a  constant  feedback 
gain  at  all  frequancies  is  not  realizabla.  Thus  tha 
bandwidth  of  the  closed  loop  system  must  still  be 
checked  to  ensure  that  the  bandwidth  is  net  excessive 
and  that  the  loop  is  rolled  off  sufficiently  where 
unmodelled  dynamics  might  affect  the  closed  loop  sta¬ 
bility  of  the  system.  The  following  robustness 
theorem  checks  this  in  a  more  formal  way. 

Robustness  Theorem  [5,6) : 

Let  G(s)  be  the  transfer  function  matrix  of  a 
proper  finite  dimensional  linear  system.  Let  G(s)  be 
its  nominal  (design)  value,  also  proper  and  let 

AG. (s)  and  JG  (s) ,  defined  by  the  relations, 

A  M 

G(s)  -  G(s)  +  AG  (s)  (49) 

A 

GCs)  -  (I+AG  (s) )G(s)  (50) 

w 

be  such  that  G(s)  and  G(s)  have  the  same  number  of 
open-loop  poles  in  the  closed  right  half  plane.  If 
the  nominal  feedback  system 

G(s)  [I+Gts))'1  -  (I+C_1(s))'1  (51) 

is  stable,  then  the  perturbed  system 

G(s)  [I+Gts)]"1  (52) 

is  stable  for  all  perturbations  satisfying 

0(I+G(jw))>  0(AGa(Ju))  VW  (53) 

or 

oCI+g'1  (jui) )  >  0  (AG  ( jw) )  VW  .  (54) 

Alternatively,  the  conditions  (53)  and  (54)  may  be 
replaced  by  (55)  and  (56)  respectively  given  as 

|  |(I+G(jw)f  l|fl>  ||AGa(}U>||  Vti)  (55) 

I!(l+G~1(ju))>‘1|r1>  l|AGM<ja>)||  VW  (56) 


Therefore  if  at  seme  frequency  uQ,  o(l+G(jw))  is  wh,re  ||.||  denotes  the  1  or  *  matrix  norms. o 

small,  there  exists  a  small  perturbation  in  the 


x 
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Computationally ,  |  |  *  |  |j  or  |  |  •  |  [^  are  cheaper  and 

easier  to  compute  than  0(*)  and  thus  are  used  in  the 
subsequent  discussion  in  place  of  the  singular  values [7J 
To  test  the  robustness  of  the  modified  LQ 
controller  for  the  multi terminal  dc  system,  several 
plots  of  the  norms  specified  in  conditions  (55)  and 
(56)  have  been  computed  with  respect  to  different  per¬ 
turbed  models.  The  nominal  model  consists  of  the 
open- loop  system  regulated  by  the  LQ  major  controller 
modified, as  previously  described,  by  passing  the  con¬ 
trol  inputs  through  washouts  with  transfer  function 

~3-^—  to  eliminate  response  to  constant  frequency 

errors.  The  perturbed  model  is  the  same  as  the  nom¬ 
inal  except  that  it  includes  additional  filtering  of 
the  control  inputs  by  an  approximately  unity  gain 
double  pole  filter  at  s=25  to  simulate  the  additional 
phase  lag  of  torsional  filters  in  the  real  system. 

The  plot  for  the  multiplicative  perturbation,  Figure 
11,  shows  that  the  perturbed  system  is  stable  and  that 
the  nominal  system  has  a  degree  of  robustness  even 
beyond  the  perturbation.  It  is  interesting  to  note 
that  the  additive  perturbation,  Figure  10,  for  the 
same  perturbed  model  exceeds  the  nominal  at  the  system 
resonances.  This  points  out  the  fact  that  the  con¬ 
ditions  (55)  and  (56)  are  only  sufficient  but  not 
necessary.  Here  in  lies  one  of  the  major  weaknesses 
of  the  singular  value  approach;  it  is  sometimes 
unnecessarily  conservative.  In  Figures  12  and  13  the 
LQ  Major  controller  is  checked  for  robustness  with 
respect  to  a  perturbed  system  that  includes  ar.  ad¬ 
ditional  machine,  with  relatively  small  inertia,  that 
is  weakly  coupled  to  the  rest  of  the  system  to  model 
intra-area  dynamics.  The  resulting  changes  in  the 
closed— loop  system  eigenvalues  are  insignificant  even 
though  the  sufficient  conditions  indicate  there  does 
exist  a  perturbation  of  the  same  size  as  the  addition¬ 
al  intra-area  dynamics  that  will  destabilize  the 
system.  Notice  also  in  the  Figures  10  and  .11  that 

| | (I+G)  1 1 |  1  and  |  j (I+G  *)  *||  *  do  not  become  small 

at  any  frequency  which  is  indicative  of  robustness. 

To  make  the  connection  with  the  classical  SISO 
gain  and  phase  margins,  under  the  LQ  major  controller 
with  washouts ,  one  of  the  system  feedback  loops  is 
broken.  The  margins  for  this  channel  are  an  infinite 
upward  gain  margin,  a  gain  reduction  margin  of  100% 
and  a  phase  margin  of  +100  degrees.  This  also  gives 
an  indication  of  how  conservative  the  norm  plots  may 
be,  since  they  take  into  account  perturbations  in  the 
system  that  may  not  be  reasonable .  Other  robustness 
norm-plots  for  the  completely  centralized  decentral¬ 
ized  case  are  similar  in  nature.  The  centralized 

case  has  the  largest  |  |  (I+G)  ||  \  as  expected,  and 

the  decentralized  the  smallest. 

Another  method  of  evaluating  the  design  is  to 
break  feedback  loops  simulating  actuator  or  sensor 
failure  or  manual  intervention  and  check  the  stability 
of  the  resulting  system.  If  the  system  remains  stable 
under  the  breaking  of  feedback  loops  the  integrity  of 
the  system  is  said  to  be  maintained.  Even  though  a 
system  is  open-loop  stable,  this  does  not  guarantee 
the  the  closed-loop  will  remain  stable  when  feed¬ 
back  loops  are  broken.  Thus  it  is  important  to  check 
the  integrity  of  the  closed-loop  system.  For  the  LQ 
major  feedback  law  with  washouts  added,  the  system 
integrity  is  maintained  for  all  possible  combinations 
of  broken  feedback  loops.  The  same  is  true  for  the 
completely  decentralized  feedback  law  with  washouts 
also. 

Finally,  a  comparison  between  the  open-and  closed- 
loop  bandwidths  is  made.  In  Figure  14  ||G(ju>)||,  the 

maximum  open-loop  frequency  response,  and 

|  |  (I+G"1(ju)))'1|  |,  the 
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maximum  closed— loop  response,  are  plotted  against  fre¬ 
quency.  This  shows  that  the  control  law  just  damps 
out  the  resonances  in  the  open— loop  system  without  af¬ 
fecting  the  rest  of  frequency  response.  Control 
energy  is  not  wasted  in  modifying  the  open— loop  system 
more  than  necessary.  This  is  related  to  the  perfor¬ 
mance  of  the  control  system  which  is  discussed  next. 

V.  CONTROL  SYSTEM  PERFORMANCE 

Evaluation  of  the  closed— loop  system  performance 
is  done  by  examining  the  system  transient  responses 
and  feedback  control  inputs  to  the  system  generated 
by  an  initial  condition  corresponding  to  the  state 
vector  for  the  open— loop  system  after  the  clearing  of 
a  three-phase  ground  fault  at  area  4.  Figure  15 
shows  the  angle  of  area  4  using  area  3  as  a  reference 
(i.e.  6^-6^)  f°r  the  open— loop,  the  LQ  centralized, 

LQ  major  and  LQ  decentralized  designs  with  the 
washouts  in  place.  Figure  16  shows  the  corresponding 
frequency  to4  for  the  same  initial  condition  and  con¬ 
trollers  . 

As  can  be  seen  from  these  figures,  all  the  con¬ 
trollers  give  approximately  the  same  performance  in 
the  transient  responses.  They  significantly  reduce 
the  settling  time  of  the  inter-area  oscillations. 

These  responses  are  typical  of  all  the  angle  and  fre¬ 
quency  transient  responses.  However,  the  control 
inputs  for  the  various  controllers  are  distinguishable. 
Figure  17  displays  the  control  signal  u2  for  the  three 

different  controllers,  again  with  the  washouts  in 
place.  As  can  be  seen,  the  centralized  feedback  gain 
uses  the  least  control  energy  indicating  a  coordinated 
control  effort.  The  decentralized  controller,  as 
expected,  use  the  most  since  it  lacks  information  to 
yield  a  coordinated  effort.  The  effect  of  the  com¬ 
munication  of  the  frequency  of  area  one  to  all  other 
areas  can  be  assessed  in  terms  of  the  difference  in 
the  control  energy  for  the  IQ  major  and  LQ  decentral¬ 
ized  cases. 


VI,  SUMMARY 

In  this  paper  some  of  the  fundamental  concepts 
of  linear  system  theory  are  exploited  to  gain  a  better 
understanding  of  the  internal  dynamics  of  the  system 
and  how  it  is  related  to  the  external  environment. 

It  is  shown  how  the  right  and  left  eigenvectors  help 
physically  interpret  the  response  of  the  system  at 
different  outputs  when  excited  from  different  inputs 
making  use  of  the  concepts  of  modal  controllability 
and  observability.  A  method  for  selecting  the  state 
weighting  matrix  to  relocate  only  specified  modes  is 
given.  The  concept  of  gain  observability  is  utilized 
in  comparing  different  feedback  structures,  with 
limited  state  feedback,  with  the  centralized  full 
state  feedback  controller.  It  is  seen  through  simu¬ 
lation  that  the  cost  of  a  decentralized  feedback 
structure  is  that  more  control  energy  is  required  to 
damp  out  oscillations  due  to  the  loss  of  communica¬ 
tion  of  important  information.  The  controllers  are 
shown  to  be  robust  with  respect  to  model  uncertainty 
by  use  of  multivariable  and  classical  robustness 
measures.  The  integrity  of  the  system  is  also  main¬ 
tained  in  the  face  of  all  possible  actuator  or 
sensor  failures.  These  considerations,  that  the 
practicing  control  engineer  must  take  into  account  in 
the  design  of  a  practical  control  system,  are  illus¬ 
trated  in  this  design  and  its  evaluation.  It  is 
important  to  emphasize  that  the  multiterminal  dc  con¬ 
troller  design  is  based  on  a  linearized  model  and 
does  not  consider  large  signal,  nonlinear  phenomena 
such  as  converter  mode  switching.  These  are  topics 
for  future  research. 
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Figure  1;  A  Power  System  with  5  DC  Terminals  and 

5  AC  Nodes,  each  representing  an  Area  or 
Group  of  Coherent  Generators. 


Figure  3:  Open-Loop  Eigenvalues  with  Nonnegative 
Imaginary  Parts. 
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Figure  4;  Ax g and  Diagrams  of  Oscillatory  Modes. 
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Figure  6:  Modal  Observability  and  Controllability 
of  a  Specified  Mode  from  Different 
Inputs  or  Outputs. 
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Figure  5:  Modal  Observability  and  Controllability 
of  Different  Modes  from  a  Specified 
Input  or  Output. 


-06  -0.7  -06  -05  -04  -03  -Q2  -0.1 


Figure  7 ;  Closed-Loop  Eigenvalues  for  Different 
Controllers. 
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Figure  10:  Robustness  for  Additive  Perturbation 
of  Additional  Torsional  Filters. 
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gure  13:  Robustness  for  Multiplicative 

Perturbation  of  Additional  Intra-Area 
Dynamics. 


gura  14 ;  Maximum  Open-Loop/Closed-Loop 

Multivariable  Frequency  Response. 


Figure  16 :  Comparison  of  Response  for 
Different  Controllers. 


Figure  17;  Comparison  of  u^  Input  for  Different 
Controllers . 
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Figure  15:  Comparison  of  Response  for 


Different  Controllers. 


